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a^ ■ ABSTRACT 

Q , Micron-sized Stardust grains that originated in ancient stars are recovered 

^^ ' from meteorites and analysed using high-resolution mass spectrometry. The most 

widely studied type of Stardust is silicon carbide (SiC). Thousands of these grains 
have been analysed with high precision for their Si isotopic composition. Here 



;_j , we show that the distribution of the Si isotopic composition of the vast majority 

of Stardust SiC grains carry the imprints of a spread in the age-metallicity distri- 
bution of their parent stars and of a power-law increase of the relative formation 
efficiency of SiC dust with the metallicity. This result offers a solution for the 
long-standing problem of silicon in Stardust SiC grains, confirms the necessity 
of coupling chemistry and dynamics in simulations of the chemical evolution of 
our Galaxy, and constrains the modelling of dust condensation in stellar winds 
as function of the metallicity. 

Subject headings: dust, extinction — meteorites, meteors, meteoroids — Galaxy: 
abundances — stars: AGB and post- AGE 
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Introduction 



A small fraction (order of 1-100 parts per million in mass) of the matrix of primitive 
meteorites is composed of Stardust grains. These grains originated in stars, were present at 
the formation of the Solar System, and have been preserved inside meteorites in their original 
form until today. Since their discovery in the late 1980s Stardust grains have been exten- 
sively analysed and employed to constrain our understanding of nucleosynthesis, mixing, 
and dust formation in stars and supernovae. Galactic chemical evolution (GCE), processing 
of dust in the interstellar r nedium, the formation of the Solar System, arid the evolution 
of meteorite parent bodies (JBernatowicz fc ZinnerlllQQTt IClayton fc Nittlerl l2004j ) . A large 
variety of minerals have been discovered as Stardust, from diamond, to silicate and AI2O3 
grains. Among them, silicon carbide (SiC) grains have been the most widely studied both 
due to their relatively large size (up to several /xm), as compared to other types of Stardust, 
and to the relatively easier separation procedure. The vast majority of the Stardust SiC 
grains recovered from meteorites show the clear si gnature of an origin in the mass-losing 
envelopes of asymptotic giant branch (AGB) stars ( IGallino et al.lll990l : iLugaro et al.l 120031 ) 
that become C rich (C>0) due to dredge-up of material from the deep C-rich layers into 
the convective envelope of the star. In C-rich conditions some C is free from the strong CO 
molecular bond and can react with Si to form SiC. Stardust SiC grains of C-rich AGB origin 
are the "mainstream" population, which comprises >93% of all Stardust SiC, and the minor 
Y and Z populations, which comprise ~1% each of all Stardust SiC. The Si compositions 
of mainstream, Y, and Z grains carry the signature of both the initial composition of their 
parent star, which is determined by GC E, and the neutro n-capture and mixing processes 
that occurred in the AGB parent stars ( iZinner et al.ll2006l ). The typical Si isotopic ratios 
of the Y and Z grains point to AGB parent stars of average metallicity c^ 1/2 and ~ 1/3, 
respectively, of solar, while the Si isotopic ratios of rnainstream grains suggest a close-to-sola r 
metallicity for their parent stars (JHoppe et al.lll997l : lAmari et al.ll200ll : IZinner et al.l 120061 ). 



A large amount of high-precision Si isotope data from Stardust SiC has been collected 
in the past 25 years (Figured]), however, their distribution is not understood. A particularly 
irksome problem is that most mainstream grains have ^^Si/^^Si and ^°Si/^®Si larger than solar 
(up to +20%) while according to GCE models t he ^^Si/^^Si and ^°Si/^^Si ratios increase with 
metallicity, which in turn increases with time (JTimmes fc Clayton! Il996[ ). The grains must 
have formed in stars of metallicity higher than solar, however, their parent stars must have 
died before the Solar System formed. Several possible explanations have been proposed 
for this app arent paradox, from a simple model of stellar migration from the inner part of 



the Galaxv (IClavton 



19971) ■ which has difficulties in reproducing the observe d distribution 



(iNittler fc Alexanderl 119991 ) , to inhomogeneities in the interstellar medium (JLugaro et al. 
19991 ). which is at odds with the correlation between the Si and Ti isotopic composition of 
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the grai ns (lNittleiil2005l) . to a starburst triggered by the merging of our Galaxy with another 
galaxy (IClaytonll2003[ ). To address this problem, we use the measured 6^^Si and 5'^°Si of 
SiC grains from AGB stars (Figure [T]) to derive the age-metallicity relation (AMR) of their 
parent stars and compare it to that observed for stars in the solar neighbourhood. 



Method 



We selected from the Presolar Grain Database (JHynes fc Gyngardll2009l ) the 2,732 main- 
stream, 133 Y, and 92 Z grains with la error bar lower than 15%o. For each grain we applied 
the following steps: 



1. We inferred the metallicity [Fe/H] (defined as the logarithm of the Fe/H ratio with 
respect to solar) of the parent star of each SiC grain from the relationship between 
S'^^Si and [Fe/H] predicted by different GCE models (Figure |2]). All the models have 
been renormalised so that at the time of the formation of the Sun, set to 8.5 Gyr, 
[Fe/H]=0 and 5^^Si=0, by definition. Before normalisation, all the models produce 
6'^^Si between —600 and —400 at [Fe/H]=0, a long-standing problem probably related 
to the rates of the nuclear reactions that produce ^^Si in core-collapse supernovae 
jHoppe et allboooh . 



2. We estimated the change in ^^''Si resulting from AGB nucleosynthesis as the distance 
A^°Si between the measured 6^^Si and the value obtained from the GCE line. We 
assumed that the best fit to the Si isotopic ratios in the Si three-isotope plot (the 
"mainstream line") shifted by —15 in 6^^Si represents the GCE of the Si isotopic ratios 
(Figured]). This shift is in agreement with the Si composition of Stardust silicate grains 
from AGB stars, which represents t he Si composition of 0-rich AGB stars as unalt ered 
by nucleosynthesis and dredge-up (JMostefaoui fc Hoppell2004j : iNguyen et al.ll2010[ ). 



We derived the parent star age from its mass as obtained from A^'^Si using the set of 
F RANEC C- rich AGB models with the neutron-c a pture cross sections of the Si isotopes 
by iGuber et al.l (120031 ) presented in IZinner et al.l ( 120061 ) , to which we added the age of 
the Sun. In these models the mass-loss rate was included using the parametrization 
given by iReimerd ( 1l975l ). For the 1.5 M0 and 2 M© models results were presented 
for different values of the associated free parameter r]=0.1, 0.3, 0.5 and we choose to 
average the results for each mass. A range of A^°Si is allowed during the C-rich phase of 
each model since A^°Si increases with the number of dredge-up episodes. This results 
in different possible masses associated to the same A^°Si, in particular when its value 
is small. We chose to remove this degeneracy by taking as the best representative of 
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each model the Si composition reported after the last computed dredge-up episode. 
Most SiC grains formed with the composition present in the envelope after the final 
few dredge-ups since the largest fraction of the envelope mass during the C-rich phase 
is lost in these final phases. For example, considering the 3 Mq model with Z = 0.02, 
more than 2/3 of the envelope mass during the C-rich phase is lost after the third-last 
TDU episode, and ro ughly 1/2 is lost after the very last TDU episode (see Table 4 of 



Straniero et al.lll997| ). Note that we did not include in our age determination th e up 



to 1 Gyr of grain residence time in the interstellar medium (JGyngard et al.ll2009l ). 



4. Since also 6"^^ Si can be marginally affected by AGB nucleosynthesis we improved the 
estimates of age and metallicity by repeating the same procedure as above with a 
new initial S'^^Simw = S'^^Si - A^^Si, with A^^Si derived from the same AGB model 
predictions used to match A^°Si. 



3. Results and discussion 

The resulting SiC AMR is plotted in Figure [3]together with the A MR derived for stars i n 



the solar neighbourhood from the Geneva-Copenhagen (G-C) survey (JHolmberg et al.ll2007l ). 



The SiC ages are affected by several uncertainties. First, there are random uncertainties 
related to the measurement errors. We made the conservative choice to plot the lower limits 
of the ages derived from adding the experimental 2a error bar to (5^°Si. Second, there are 
random errors related to the possible effect of inhomogeneities in the i nterstellar medium 



which may ch ange the Si composition of any given parent stars by ~ 50%c (JLugaro et al.lll999 



Nittlerll2005l ). These are not possible to be evaluated. Assuming that the silicon isotopic 
distribution is relatively symmetric, the age and the metallicity calculated for grains at the 
peak of the distribution should be reliable. Third, there are systematic uncertainties related 
to the choice of the line taken to represent the GCE in the Si isotope plot (Figure [T]) and 
to the AGB model predictions. The uncertainty related to the GCE line would most likely 
result in smaller stellar ages as the line could be shifted further away from the mainstream 
line than what we have assumed and still be in agreement with the silicate data. Instead, it 
is not known if the uncertainty related to the AGB models would result in smaller or larger 
stellar ages since both larger and smaller A^°Si for a given stellar mass are possible within 
the uncertainties wrought by, e.g., the mass- loss rate, the efficiency of the dredge- up of the 
deep layers of the stars into the convective envelope, and the neutron-capture cross section 
of the Si isotopes. 

Nonwithstanding the uncertanties discussed above, relatively, the ages derived for the 
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parent stars of the Y grains are similar to those of the mainstream grains, while those derived 
for the Z grains cover a much narrower range indicating that Z grains have on average parent 
stars of higher mass than mainstream and Y grains. This result supports proton captures 
at the base of the convective envelope (also known as "hot bottom burning" ) as the process 
responsible for lowering the ^^C/^^C ratios in the Z grains (ranging from 20 to 100) with 
respect to those observed in the Y grains (> 100, by definition). If this interpretation 
is correct, GCE rn odels are not required to match the 6'^^Si versus [Fe/H] relationship of 
Zinner et al.l ( 120061 ) (Figure [2]), which was derived using AGB models of fixed mass lower 
than ~ 3 M0. On the other hand, one may wonder why there are no SiC grains from low- 
metallicity and low-mass parent stars. A possibility is that the very high C/0 ratio reached 
in th ese stars (up to 2 - 30) may favour production of amorphous carbon dust rather than 



SiC fjSloan et al 



200 



The SiC [Fe/H] distribution is determined by the steepness of the S'^^Si versus [Fe/H] 
relationship, which varies with the choic e of the core-collapse su pernova (SNII) yields and 
the GCE model (Figure [2]). The GEtool ( Fenner fc GibsonI l2003l ) chemical evolution model 
was computed with dual infall (i.e., a rapid formatio n of the halo, followe d by subsequent, 
protracted, disk formation), initial mass function from lKroupa et al.l (119931 ). and a Schmidt- 
Kennicutt star formation prescription, and it is tuned to recover the gas and stellar abun- 
dances and radial surface densities in the Milky Way. In Figure [3] we plot two SiC AMRs 
obtained using the most and the least steep (5^^Si versus [Fe /H] relationships from the 



GEtool simulations. 



and iKobavashi et al.l ( 



^'hese corresponds to using the yields by IWoosley fc Weaver! (119951 ) 
2006|), resp ectively. The GEtool mo del computed using the SNII 



yields from iKobayashi et al.l (120061 ) and that computed using IChieffi fc Limongil (120041 ) pro- 
duce a very similar SiC AMR with a [F e /H] sp read of a factor of three, wit hin that observed 
in the G-C survey. iTimmes fc Claytoru (119961 ) and IKobayashi et al.l (120 111 ) obtained a much 
steeper and much flatter, respectively, 5^^ Si and [Fe/H] relationship than the GEtool models 
presented here. The difference depends on many of the ingredients of the GCE simulation, 
e.g., the initial mass function and the infall scheme, which also affect the theoretical AMR 
predicted directly by the GCE models. Interestingly, GCE models that predict a steeper 
5^^Si versus [Fe/H] relationship also predict a flatter AMR. As a consistency check, in the 
right panel of Figure [31 the AMRs predicted by the different GCE models considered here are 
compared to the G-C stars. None of the models can recover the observed AMR. A better 
match may be found by recent models that corisider dynamics together wit h the chemical 
evolution in the Galaxy (IKobayashi fc Nakasatdl201ll : iPilkington et al.ll2012n . However, we 
do not consider them here as they have not been extended yet to include the evolution of 
isotopic abundances. Overall, the SiC grain data confirm the result of the G-C survey that 
stars exist with ages older than the Sun and metallicities higher than the Sun. 
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To make a quantitative comparison of the metallicity distribution function (MDF) de- 
rived for the SiC grain parent stars and for the G-C stars we restrict ourselves to the main- 
stream SiC grains as they represent the least biased sample. We removed the Y and Z grains 
because their numbers are probably overestimated due to specific searches dedicated to iden- 
tifying these types of grains. In the left panel of Figure H] we compare the MDF obtained 
from the SiC AMR to that obtained from the G-C survey. Because the la error bar on the 
grain [Fe/H] that derives from the experimental uncertainty of S'^^Si is much smaller (<0.03 
dex) than that of the stellar [Fe/H] (~0.1 dex), for better comparison we convolved the grain 
[Fe/H] data with a Gaussian of a=0.1. The main difference between the grain and the stellar 
MDF is that the mean metallicity of the parent stars of the mainstream SiC grains is ~ 50% 
higher than that of the G-C stars. We interpret this as a selection effect between the grain 
and the star samples indicating that formation of SiC dust is favored by higher metallicity. 
We define the SiC relative formation efficiency (RFE) as the ratio between the normalised 
number of mainstream SiC grains and of G-C stars in each metallicity bin and plot these 
values in the right panel of Figure HI The SiC RFE is unitless and its values are not absolute, 
but have a meaning only when considered relatively to each other. We infer that the SiC 
RFE can be represented by a power law in metallicity. The relationship is well defined only 
for values of [Fe/H] between ±0.3 dex because in this range there are sufficient numbers of 
both mainstream SiC grains and G-C stars to make their ratio statistically meaningful. 

An important systematic uncertainty in the derivation of the SiC RFE is related to the 
renormalised value of 5^^Si=0 at [Fe/H]=0 (Item 1 of Sec. |2]) Such renormalisation may easily 
be a few percent wrong, in particular due to the effect of inhomogeneities in the interstellar 
mediur n, which prod uce variations in (5^^Si of the order of 70%o for the same [Fe/H] within 



±0.01 ( jNittlerl |2005| ). If 5^^Si is renormalised to a positive value instead of zero, clearly 
the grain MDF becomes closer to that of the G-C survey and the derived SiC RFE is less 
steep. In Figure S] we also present the results obtained by setting 5^^Si=50%o at [Fe/H]=0. 
The SiC RFE still increases with the metallicity, though the increase is less pronounced and 
disappears between 0<[Fe/H]<0.1. We consider this example as an upper limit: a choice 
of 5^^Si~30%o at [Fe/H]=0 is probably more realistic being the valu e shown by the largest 



number of SiC, according to Figure 13 of iNittler fc Alexander! ( 120031 ) 



4. Conclusions 

The results presented here confirm that the relationship between age and metallicity in 
the Galaxy is relatively flat and that a spread of metallicities is present for each age. This can- 
not be recovered by traditional GCE models and requires a more sophisticated approach in- 
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cluding coupling of Ga lactic dynamics and chemical evolution ( iKobayashi &: Nakasatdl201ll : 
Pilkington et al.l |2012| ). Scattering and radial migration must have pl ayed an i mpor tant 



role in determining the properties of stars in the solar neighbourhood (IClavton 



supported by recent observational studies (e.g. iBoeche et al. II2013I : iRamirez et al 



997) 



as 



2013|) 



We have interpreted the shift of the SiC MDF to higher metallicities than the G-C survey 
as a selection effect and derived that the SiC relative formation efficiency increases with 
the stellar metallicity as a power law. This result is in qualitative agreement with Spitzer 



obser vations of C-rich AGB stars in the Large and Small Magellanic Clouds (jSloan et al. 



20081 ). which indicate that the mid-infrared emission from SiC and silicate dust decreases with 
the metallicity, while the emission from amorphous carbon dust does not. On the basis of 
this evidence, we tentativ ely predict that the MDF of Stardust silicates ( iMostefaoui fc Hoppe 
2004 : iNguyen et al.ll2010l ) should be simil ar to that of SiC gra ins, while the MDF of Stardust 
graphite grains from C- r ich A GB stars (jJadhav et al.l |2008| ) should be similar to that of 
the G-C survey. iNittlerl ( l2009l ) conducted a similar exercise to ours based on the O ratios 
in Stardust oxide grains and found evidence for the existence of a moderate AMR. If the 
i8q^i6q ti-q;^[q is g^ good indicator of the stellar metallicity, we expect a correlation between 
i8q^i6q g^-^^ S'^^Si, which is not shown by the Stardust silicate gra ins, but could be masked 
by the effect of dilution with normal material ( INguyen et al.ll2010l ). It should also be noted 
that Stardu st oxide and silic ate grains are more likely to come from lower-mass stars than 
SiC grains ( IGail et al.ll2009l ). which may result in a different AMR. It will be possible to 
statistically investigate these issues in the future when more high-precision Si data for these 
types of Stardust are available. 
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Fig. 1. — Distribution of the Si isotopic compositions of SiC mainstream, Y, and Z grains 
with la error bar lower than 15%o represented using the (5^^'^°Si notation, i.e., the permil 
variation of ^^'^°Si/^^Si with respect to the solar ratios. The solid line is the best fit through 
the mainstream SiC data ("mainstream line") of slope 1.31 and intercept — 15.9%o. The 
dashed line is the mainstream line shifted by — 15%o in 5^°Si. This is taken to represent the 
GCE of the Si isotopic ratios, with 5^^'^° Si increasing with the metallicity. 
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Fig. 2. — T he relationship between ^ ^^Si and [Fe/H] derived 
models: TC96 (ITimnies fc Claytoru Il996l ). KKUll (JKobayashi et al. 



Tom different GCE 



,, 120111). and GEtool 

jFenner fc Gibson! boosh using different S NII yields: W W95 JWooslev fc Weaveil llQQsh. 



CL04 (IChieffi fc Limongil I2004J ) , and K06 (jKobayashi et al.l 120061 ) ; or by simply assuming 
that the abundances of ^^Si and ^''Si scale with the metallicity, while that of ^*Si is a en- 
hanc ed such that ^^ Si i s 1/8 of its solar abundance when the metallicity is 1/10 of solar: 
Z06 (IZinner et al.l 120061) .This c hoice produces a (5^^Si and [Fe/H] relationship close to that 
obtained by IZinner et al.l (120061 ) comparing the composition of Z grains to AGB models of 
fixed mass and variable metallicity. 
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Fig. 3. — The age-metallicity relation (AMR) derived for the parent stars of SiC mainstream, 



Y, an d Z grains using GEtoo l with two different sets of SNII yields: K06 (JKobayashi et al. 



20061 . left panel) and WW95 ( IWoosley fc Weaverlll995l . middle panel). Note that the plotted 
ages are lower limits based on the experimental 2a error bar. The upper limits are unde- 
termined (because of negative A^°Si) and >13 Gyr for 54% and 26% of the mainstream 
grains, respectively, while they are <6 Gyr for 86% of the Z grains. Selected initial stellar 
masses are indicated at the top of the panels in correspondance to their ages. Note that 
only AGB stars in the mass range between ^1.5 Mf^ and ~4 M^r^ are ex pected to become C 
rich and produce SiC grains (JGroenewegen et al.lllQQSuGail et al.ll2009l ). The corresponding 
age limits are highlighted by the two vertical thin black lines, however, we did not remove 
the points outside this range because (i) they are lower limits and (ii) they could shift in- 
side the allowed range when considering the several uncertainties discussed in the text. The 
right panel shows the AMR obtained from the Geneva-Copenhagen (G-C) survey for 2037 
stars in the solar neighbourhood with age uncertainties lower than 25%, as compared to the 
AMRs predicted by the different GCE mod els: GEtool (solid line, w hich is independent of 
the choice of the yield s) ; TC96 (dashed line. iTimmes fc Claytonlll996l ): KKUll (dotted line, 
Kobayashi et al.ll201ll ). 
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Fig. 4. — The metallicity distribution function (MDF) of the parent star s of mainstream SiC 
grains derived u sing GEtool and two dif ferent sets of SNII yields: WW95 fjWoosley fc Weaver 
19951 ) and K06 (JKobayashi et al.l 120061 ) . is compared to that obtained for stars in the solar 
neighbourhood from the Geneva- Copenhagen (G-C) survey (left panel). The SiC relative 
formation efficiency (RFE) of as function of the [Fe/H] is shown in the right panel. The 
thin lines present the results obtained renormalising (5^^Si=50%o at [Fe/H]=0 and the black 
dashed line is a proposed power-law fit. 



